********* Genera la Base de datos Modified_BaseContrasteInteres.dta *******

*********************************************
****Calculos con base en ENLACE 2012-2013****
*********************************************

clear
set more off

use "$basein/Datos adicionales/Enlace_2013.dta", replace

***Traigo base de tratamientos***

merge m:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

***Dejo solo las observaciones que necesito***

keep if _merge==2 | _merge==3

******************
***Base estrato***
******************
gen school_id=cct

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

keep if _estratos==3

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"



***Genero cct numerico***

egen cct_num=group(cct)

***Regresiones simples***

gen PEC_Plus=0
replace PEC_Plus=1 if treatgroup2=="PEC Plus"

gen PEC=0
replace PEC=1 if treatgroup2=="PEC"

***Estandarizo variables resultado***

foreach x in 3 4 5 6 {

egen esp_std_`x'=std(cal_esp) if grado=="`x'"
egen mat_std_`x'=std(cal_mat) if grado=="`x'"

}

gen esp_std=.
gen mat_std=.

foreach x in 3 4 5 6 {

replace esp_std=esp_std_`x' if grado=="`x'"
replace mat_std=mat_std_`x' if grado=="`x'"

}

***Medias simples***

table grado treatgroup2, c(mean mat_std)
table grado treatgroup2, c(mean esp_std)

save "$base_out/Enlace_2013ftTratamientos.dta", replace

*********************************
****Calculos con base en 911*****
*********************************

clear
set more off

use "$basein/Datos adicionales/911_interes.dta", replace

***Traigo base con variables de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

***Dejo solo las variables de interes***

keep if _merge==3

***Traigo identificador de estrato para tablas***

rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Regresiones simples***

gen PEC_Plus=0
replace PEC_Plus=1 if treatgroup2=="PEC Plus"

gen PEC=0
replace PEC=1 if treatgroup2=="PEC"


***Estad'isticas simples***

table treatgroup2, c(mean f_T_estu mean f_T_profe mean f_PT_ratio)
table treatgroup2 if Nivel=="Primaria", c(mean f_T_estu mean f_T_profe mean f_PT_ratio)

foreach var in f_T_estu f_T_profe f_PT_ratio{
	sum `var'
	matrix tempm=r(mean)
	local tempm=string(tempm[1,1], "%9.2fc")
	file open newfile using "$tables/`var'_mean.tex", write replace
	file write newfile "`tempm'"
	file close newfile
}
		
*******************************************************
****Cálculos marginacion ruralidad e infraestructura***
*******************************************************


clear
set more off
use "$basein/Datos adicionales/base_Contraste.dta"

drop _merge

rename school_id cct

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

***Traigo cambio de director***


rename cct school_id

merge 1:1 school_id using "$basein/Datos adicionales/Change_Director_2015_2018.dta", gen(_cambio)

***Traigo identificador de estrato para tablas***

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"


***Regresiones simples***

gen PEC_Plus=0
replace PEC_Plus=1 if treatgroup2=="PEC Plus"

gen PEC=0
replace PEC=1 if treatgroup2=="PEC"

***Genero variable de marginación***

gen marginales=.
replace marginales=0 if marginacion!=""
replace marginales=1 if marginacion=="Muy alta o Alta"

***Estad'isticas simples***

table treatgroup2, c( mean marginales mean urbana mean infra mean cambio)

foreach var in marginales urbana cambio{
	sum `var'
	matrix tempm=r(mean)*100
	local tempm=string(tempm[1,1], "%9.2fc")
	file open newfile using "$tables/`var'_mean.tex", write replace
	file write newfile "`tempm'"
	file close newfile
}

save "$base_out/base_ContrasteftChangeDirector.dta", replace
	
***********************
**Niveles Planea 2015**
***********************

clear
set more off
use "$basein/Planea/Resultados_planea_2018_nivelescuela.dta"

drop in 1/2

drop if entidad=="CHIAPAS"

rename clavedelaescuela school_id

merge 1:1 school_id using "$basein/Datos adicionales/base_Contraste.dta", gen(planea)

keep if planea==2 | planea==3

rename school_id cct
drop _merge

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

**Genero tratamiento**

gen PEC_Plus=0
replace PEC_Plus=1 if treatgroup2=="PEC Plus"

gen PEC=0
replace PEC=1 if treatgroup2=="PEC"

gen Control=0
replace Control=1 if treatgroup2=="Control"

**Genero variables para la tabla**

gen Mat_low50=mati+matii
gen Mat_Hig50=matiii+mativ
gen Leg_low50=lengi+lengii
gen Leg_Hig50=lengiii+lengiv

save "$base_out/ModifiedResultados_planea_2018_nivelescuela.dta",replace

*********************************************
**Graficos relación entre DWMS y resultados**
*********************************************

clear
set more off
use "$basein/Planea/Resultados_planea_2018_nivelescuela.dta"

drop in 1/2

drop if entidad=="CHIAPAS"

rename clavedelaescuela school_id

merge 1:1 school_id using "$basein/Datos adicionales/base_Contraste.dta", gen(planea)

keep if planea==2 | planea==3

rename school_id cct
drop _merge

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

**Genero tratamiento**

gen PEC_Plus=0
replace PEC_Plus=1 if treatgroup2=="PEC Plus"

gen PEC=0
replace PEC=1 if treatgroup2=="PEC"

gen Control=0
replace Control=1 if treatgroup2=="Control"

**Genero variables para la tabla**

gen Mat_low50=mati+matii
gen Mat_Hig50=matiii+mativ
gen Leg_low50=lengi+lengii
gen Leg_Hig50=lengiii+lengiv

**Traigo base de baseline DWMS**

merge 1:1 cct using "$basein/Datos adicionales/DWMS_fase1.dta", gen(DWMS)

keep if DWMS==3

**Hago graficas

twoway (scatter Mat_Hig50 index_nuevo) (lfit Mat_Hig50 index_nuevo, lwidth(thick) ), graphregion(color(white)) xtitle("Management index") ytitle("Share of students in top 50% Planea (Math)") legend(off)

twoway (scatter Leg_Hig50 index_nuevo) (lfit Leg_Hig50 index_nuevo, lwidth(thick) ), graphregion(color(white)) xtitle("Management index") ytitle("Share of students in top 50% Planea (Lang)") legend(off)

************************
****Datos reescuchas****
************************

clear
set more off

use "$basein/Datos adicionales/base_Contraste.dta"

drop _merge

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

drop _merge

merge 1:1 cct using "$basein/Datos adicionales/Datos_3_Reescucha.dta"

merge 1:1 cct using  "$basein/Datos adicionales/911_interes.dta", gen(_911)

keep if _911==3


***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

***Traigo cambio de director***


merge 1:1 school_id using "$basein/Datos adicionales/Change_Director_2015_2018.dta", gen(_cambio)

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if _merge==3

***Genero variable de marginación***

gen marginales=.
replace marginales=0 if marginacion!=""
replace marginales=1 if marginacion=="Muy alta o Alta"


***Estadísticas simples***
tabstat urbana marginales infra f_T_estu f_T_profe f_PT_ratio cambio, by(Reescuchado)

save "$base_out/base_ContrasteDatosReescuchas.dta", replace


***************************************
*****Ahora reescuchas en planea 2018?***   *decia2015*
***************************************

clear
set more off
use "$basein/Planea/Resultados_planea_2018_nivelescuela.dta"

drop in 1/2

drop if entidad=="CHIAPAS"

rename clavedelaescuela school_id

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

**Traigo base de baseline DWMS**

merge 1:1 cct using "$basein/Datos adicionales/Datos_3_Reescucha.dta", gen(DWMS)

***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if DWMS==3

**Genero variables para la tabla**

gen Mat_low50=mati+matii
gen Mat_Hig50=matiii+mativ
gen Leg_low50=lengi+lengii
gen Leg_Hig50=lengiii+lengiv

***Estadísticas simples***
tabstat Mat_low50 Leg_low50, by(Reescuchado)

save "$base_out/Resultados_planea_2018_nivelescuelaReescuchas.dta",replace
**********************************
**********************************
***Respondo o no responde valora**
**********************************
**********************************

**************************
***Observacion de clase***
**************************

clear
set more off

use "$basein/Datos adicionales/base_Contraste.dta"

drop _merge

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

drop _merge

merge 1:1 cct using "$basein/Instrumentos valora/Observacion_indices.dta"

drop if _merge==2

merge 1:1 cct using  "$basein/Datos adicionales/911_interes.dta", gen(_911)

keep if _911==3


***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Traigo cambio de director***


merge 1:1 school_id using "$basein/Datos adicionales/Change_Director_2015_2018.dta", gen(_cambio)


***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if _merge==3

***Genero variable de marginación***

gen marginales=.
replace marginales=0 if marginacion!=""
replace marginales=1 if marginacion=="Muy alta o Alta"


***Estadísticas simples***
tabstat f_T_estu f_T_profe f_PT_ratio urbana marginales infra  cambio, by(Reescuchado)

save "$base_out/base_ContrasteObsClase.dta", replace

**************************
***Herramienta de clase***
**************************

clear
set more off

use "$basein/Datos adicionales/base_Contraste.dta"

drop _merge

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

drop _merge

merge 1:1 cct using "$basein/Instrumentos valora/Herramientas_indices.dta"

drop if _merge==2

merge 1:1 cct using  "$basein/Datos adicionales/911_interes.dta", gen(_911)

keep if _911==3


***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Traigo cambio de director***


merge 1:1 school_id using "$basein/Datos adicionales/Change_Director_2015_2018.dta", gen(_cambio)


***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if _merge==3

***Genero variable de marginación***

gen marginales=.
replace marginales=0 if marginacion!=""
replace marginales=1 if marginacion=="Muy alta o Alta"


***Estadísticas simples***
tabstat f_T_estu f_T_profe f_PT_ratio urbana marginales infra  cambio, by(Reescuchado)

save "$base_out/base_ContrasteHerramientaClase.dta",replace

***********************************
***********************************
*****Ahora valora en planea 2015***
***********************************
***********************************

**************************
***Observacion de clase***
**************************

clear
set more off
use "$basein/Planea/Resultados_planea_2018_nivelescuela.dta"

drop in 1/2

drop if entidad=="CHIAPAS"

rename clavedelaescuela school_id

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

**Traigo base valora**

merge 1:1 cct using "$basein/Instrumentos valora/Observacion_indices.dta", gen(_valora)

drop if _merge==2

***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if _valora==3

**Genero variables para la tabla**

gen Mat_low50=mati+matii
gen Mat_Hig50=matiii+mativ
gen Leg_low50=lengi+lengii
gen Leg_Hig50=lengiii+lengiv

***Estadísticas simples***
tabstat Mat_low50 Leg_low50, by(Reescuchado)

save "$base_out/Resultados_planea_2018_nivelescuelaObsClase.dta",replace

**************************
***Herramientas de clase**
**************************

clear
set more off
use "$basein/Planea/Resultados_planea_2018_nivelescuela.dta"

drop in 1/2

drop if entidad=="CHIAPAS"

rename clavedelaescuela school_id

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

**Traigo base valora**

merge 1:1 cct using "$basein/Instrumentos valora/Herramientas_indices.dta", gen(_valora)

drop if _merge==2

***Traigo identificador de estrato para tablas***
rename cct school_id

merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estratos)

drop if estrato==13 | estrato==15

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Genero variable de reescucha***

gen Reescuchado=0
replace Reescuchado=1 if _valora==3

**Genero variables para la tabla**

gen Mat_low50=mati+matii
gen Mat_Hig50=matiii+mativ
gen Leg_low50=lengi+lengii
gen Leg_Hig50=lengiii+lengiv

***Estadísticas simples***
tabstat Mat_low50 Leg_low50, by(Reescuchado)

save "$base_out/Resultados_planea_2018_nivelescuelaHerrClase.dta",replace

**********************************************************************
***Tabla general participantes en el experimento vs. otras escuelas***
**********************************************************************

clear
set more off

use "$basein/Datos adicionales/base_Contraste.dta"

drop _merge

rename school_id cct

***Traigo bases de interes***

merge 1:1 cct using "$basein/Datos adicionales/Tratamientos_completos.dta"

drop _merge

merge 1:1 cct using "$basein/Datos adicionales/Datos_3_Reescucha.dta"

merge 1:1 cct using  "$basein/Datos adicionales/911_interes.dta", gen(_911)

******************
***KEEP´ONLY PEC VS PEC PLUS******
******************

drop if state==10
drop if treatgroup2=="Control"

***Estadisticas simples de escuelas***

gen evaluacion=1 if treatgroup2!=""
replace evaluacion=0 if evaluacion==.

***Identifico estados con participacion en al evaluacion***

gen estado=substr(cct,1,2)
encode estado, gen(state1)

gen participacion=0
replace participacion=1 if estado=="10" | estado=="11" | estado=="15" | estado=="17" | estado=="21"  | estado=="27" | estado=="29"

tabstat f_T_estu f_T_profe f_PT_ratio  , by( evaluacion) stats(mean) col(stats) nototal

tabstat f_T_estu f_T_profe f_PT_ratio if participacion==1 , by( evaluacion) stats(mean) col(stats) nototal


rename cct clavecct

drop id_localidad

merge 1:1 clavecct using "$basein/Datos adicionales/localidad_escuela.dta", gen(_localidad)

keep if _localidad==1 | _localidad==3

***Traigo marginacion***

merge m:1 id_localidad using "$basein/Datos adicionales/marginacion_2010.dta", gen(_marginacion)

keep if _marginacion==1 | _marginacion==3

drop urbana

***Traigo urbanidad**

merge m:1 id_localidad using "$basein/Datos adicionales/urbanas.dta", gen(_urbanas)

***Genero las variables de interés**

gen Urbanas=0
replace Urbanas=1 if _urbanas==3 

gen Alta_marginacion=0
replace Alta_marginacion=1 if gm_2010=="Alto" | gm_2010=="Muy alto"

***Genero estadisticas para la tabla**

tabstat   Urbanas, by( evaluacion) stats(mean) col(stats) nototal

tabstat   Urbanas if participacion==1 , by( evaluacion) stats(mean) col(stats) nototal

************* original values ***********
gen math25=exp(ln_math25)
gen lang25=exp(ln_lang25)
gen mathlow25=exp(ln_mathlow25)
gen langlow25=exp(ln_langlow25)
gen math50=exp(ln_math50)
gen lang50=exp(ln_lang50)
gen stu=exp(lstu)

rename clavecct school_id
merge m:1 school_id using "$basein/Datos adicionales/estratos.dta", gen(_estr)
rename school_id clavecct

merge 1:m clavecct using "$base_out/Planea2015_Escuelas.dta",gen(_Planea15)

collapse (first) marginacion (first) treatgroup2  (mean) f_T_estu (mean) f_T_profe (first) estado (first) state (mean) participacion (first) gm_2010 (mean) urbano (mean) Urbanas (mean) Alta_marginacion (mean) estrato (mean) ENT (mean) mathI (mean) mathII (mean) mathIII (mean) mathIV (mean) langI (mean) langII (mean) langIII (mean) langIV (mean) evaluacion (first) state1 (mean) f_PT_ratio , by(clavecct)

gen f_studteach_ratio=f_T_estu/f_T_profe
label variable f_studteach_ratio "# of students/# of teachers"

gen expEdoMex=0
replace expEdoMex=1 if state==9
replace expEdoMex=2 if state==10


save "$base_out/Modified_BaseContrasteInteres.dta", replace
